Cooper-pair propagation and superconducting correlations in graphene 
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We investigate the Cooper-pair propagation and the proximity effect in graphene under conditions 
in which the distance L between superconducting electrodes is much larger than the width W of the 
contacts. In the case of undoped graphene, supercurrents may exist with a spatial decay proportional 
to W^/L^. This chan ges upon doping into a behavior, opening the possibility to observe a 

supercurrent over length scales above 1 /xm at suitable doping levels. We also show that there is 
in general a crossover temperature T* ~ vp/ksL that marks the onset of the strong decay of the 
supercurrent, and that corresponds to the scale below which the Cooper pairs are not disrupted by 
thermal effects during their propagation. 



The recent interest in the physics of graphene has 
arisen from the observation of a number of novel elec- 
tronic properties P, which are the consequence of the 
relativistic-like character of the electron quasiparticles. 
This is certainly the case of the anomalous quantization 
of the Hall conductivity!^, dj IHl' ^ ^ exis- 
tence of a finite lower bound in the conductivity at the 
charge neutrality point 0, H, 0, [1| • The graphene system 
is unique in that the low-energy excitations have conical 
dispersion around discrete Fermi points, being therefore 
governed by a Dirac equation for massless chiral particles. 
The appearance of an additional pseudo-spin quantum 
number intrinsic to the Dirac spectrum has led to propose 
other unconventional effects, like the selective transmis- 
sion of electrons through n-p junctionsfo'] or a peculiar 
form of Andreev reflection at the metal-superconductor 
interface in graphene [lo|. 

Recently the proximity effect has been investigated in 
graphene in two different experiments 11 , l2| ■ In both of 
them a 2D carbon layer has been contacted with super- 
conducting electrodes, though the size and aspect ratio 
of the respective samples have been quite different in the 
two cases. Thus, while the separation of the electrodes 
in the experiment of Ref. [llj appears to be of the order 
of ~ 0.5 /xm, the distance between them in the report of 
Ref. seems to have a minimum value of about 2.5 /im. 
This may explain why a supercurrent has been measured 
between the electrodes in the first case, while the most 
relevant observation in the second experiment has been 
an abrupt drop in the resistance, below the critical tem- 
perature of the superconducting electrodes. 

This crucial dependence on the distance between the 
electrodes may be however somewhat surprising, taking 
into account that supercurrents have been measured in 
carbon nanotubes placed between superconducting con- 
tacts, with nanotube lengths as large as « 1.7 /im [13]. It 
becomes then pertinent to analyze the proximity effect in 
graphene, specially for large separation between the su- 
perconducting electrodes, in order to find out the features 
related to the 2D character of the material. The Joseph- 



son effect has been studied before in graphene strips with 
length L small relative to their width and the supercon- 
ducting coherence length|lj]. Here we will face the op- 
posite situation, in which the relevant signatures are dic- 
tated by the propagation of the Cooper pairs in graphene 
and its dependence on thermal and interaction effects. 

We investigate then the proximity effect in graphene, 
under conditions in which the distance L between the su- 
perconducting electrodes producing the Cooper pairs is 
much larger than the width W of the contacts. More 
specifically, we will consider that the hamiltonian for 
graphene with superconducting contacts along the coor- 
dinates xi = Q and X2 — L is given by 

H = VF J cPr *^'')^(r)cr(°) • d ^i^\v) 

+ *^"^^(^j,2/)*5,..(2:,,y) -|-h.c(l) 

cr(°) being different sets of Pauli matrices for a — 1,2 0]. 
In the above expression, a sum is implicit over the spin 
index a, as well as over the index a accounting for the two 
different valleys and corresponding Dirac spinors at 
opposite K points in graphene. We will consider that 
the coupling to the electron fields ^si and '^S2 in the 
superconducting electrodes takes place along a segment 
of width W in each case. For our purposes, it will be 
enough to describe the superconductors in terms of the 
order parameter A and the normal density of states N. 

We are going to deal in particular with the case in 
which the time of propagation of the Cooper pairs be- 
tween the contacts is much larger than 1/|A|. This 
implies, equivalently, that the distance L has to be 
much larger than the superconducting coherence length 
^. Then, a supercurrent may arise mainly from processes 
in which the Cooper pairs tunnel from one of the su- 
perconductors to the 2D layer, propagating to the other 
superconducting contact. The Josephson current Is can 
be computed as the derivative of the free energy with 
respect to the difference x between the respective phases 
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Xi and X2 of the order parameters in the two supercon- 
ductors, Is = —2edT/dx- Under the assumption of a 
relative large |A|, we can make the approximation 

(2) 

where the statistical averages, at temperature T, are 
taken over ordered products with respect to imaginary 
time T. Then, the maximum supercurrent (critical cur- 
rent) becomes to lowest order in the tunneling amplitude 
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In deriving Eq. ([3]) we have already assumed that L ^ 
W, averaging over the separation of the electrons in the 
pair tunneling along the superconducting contact. 

The determination of the critical current is there- 
fore reduced to the evaluation of the propagator of the 
Cooper pairs in the graphene layer. The Fourier trans- 
form of this object, that we will denote by D{]i, iZJ), can 
be obtained from the standard diagrammatics for Dirac 
fermions. Here we stress that, for the sake of preserving 
the relativistic-like invariance, it is convenient to regular- 
ize the diagrams at high energies by using a method that 
maintains the space-time symmetry of the theory, like the 
analytic continuation in the number of dimensions [l5| . In 
particular, the Cooper-pair propagator can be obtained 
in the noninteracting theory from the convolution of two 
Dirac propagators, and at T = it turns out to be 
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(4) 



The temperature dependence of this object can be ob- 
tained using the Matsubara formalism. In the static limit 
we get 
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X ln[2 cosh(wj.|k|v/a;(l -a;)/2fci3T)](5) 

The knowledge of the Cooper-pair propagator at zero 
frequency is enough to compute the critical current in Eq. 
([3]). Interaction effects can be incorporated by summing 
up the different perturbative orders obtained by iteration 
of the scattering of the particles in the pair. If we take 
an average of the potential V between the particles, we 
have that the Cooper-pair propagator can be represented 
in the T-matrix approximation by 



i:>W(k,^) 



(6) 



We remark anyhow that the interaction effects cannot af- 
fect significantly the propagation of the Cooper pairs at 
low temperatures. When there is exchange of the valleys 
in the scattering of the electrons in the pair, the Coulomb 
potential Vc becomes suppressed by a large momentum- 
transfer 2kF, down to Vc ~ /2kp. At the small values 
of |k| relevant for the long-distance regime, the denom- 
inator in Eq. ([6|) remains then very close to 1. In the 
opposite case of no valley exchange, we observe that at 
low momentum-transfer VcD'^^^ has a relative strength 
of the order of ^ /vp- At low energies, this effective 
coupling is strongly renormalized, down to a value which 
may be about one order of magnitude below the nomi- 
nal coupHng[l^. We have checked that, in practice, this 
makes very small the difference between computing the 
critical current ^ from either D(k,0) or Z)'^''^(k, 0), at 
the relevant temperatures in the experiments (of the or- 
der of - 1 K). 

It becomes convenient to factor out from ([3]) the rel- 
ative tunnel conductances of the interfaces, which are 
given each by the constant quantity Nt^W/vp [17|]. We 
will deal then with the intrinsic 2D dependence of the 
critical current, which becomes in terms of the Fourier 
transform of the Cooper-pair propagator 
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where kc is a short-distance cutoff, that we take of the 
order of ~ 1 nm~^. 

The critical current computed from ^ displays a num- 
ber of peculiar features. From the scaling of Z)^"^ (k, 0) in 
the zero-temperature limit, it is easily seen that Ic^^\o) 
must decay as ~ vpW^ /L^ at large separation between 
the contacts. Another important property is that there 
is always a crossover temperature T* ~ vp/ksL which 
marks the onset of a power-law decay of the supercurrent 
at high temperatures. This is shown in Fig. [Tl which rep- 
resents the behavior of the critical current as a function 
of T for L = 5 X 10^ /k^ and W = 10^ /k^ . The depen- 
dence of Ic'^'^^ (T) is qualitatively similar to what is found 
in carbon nanotubes with long separation between super- 
conducting contacts [3 [l^. The crossover temperature 
corresponds to the scale at which the Cooper pairs fail 
to propagate efficiently between the electrodes, as they 
become disrupted along the way by thermal effects. 

The novel behavior of the supercurrents in graphene 
with respect to that in the carbon nanotubes comes 
from the different scaling with length L at low temper- 
atures. This is given by a power-law behavior in 
graphene, instead of the 1/L scaling of a noninteract- 
ing ID system^lSi]. We may consider for instance the 
sample described in Ref. with an approximate sepa- 
ration between the electrodes L w 2.5 /xm. The plot in 
Fig. [T] gives the theoretical values of the critical current. 
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FIG. 1: Plot of the critical current /f"'(r) (in units of 
lO-^evpkc ~ 1.2 fj,A) for L = 5 X lO^kc (~ 2.5 fim) and 
W = lOVfcc (~ 50 nm). 



in units of 1Q~^ evpkc « 1.2 /lA. We observe that, in 
the low-temperature regime to the left of the crossover, 
the expected critical currents should have a magnitude 
Ic'^'^^ ^ 10^^ nA. Such an small scale may explain the 
difficulty to establish a supercurrent in graphene when 
there is a separation of the order of microns between the 
superconducting electrodes. 

On the other hand, the crossover shown in Fig. [T]has a 
remarkable correspondence with the abrupt drop in the 
resistance reported in Ref. 



12. This feature has been 



observed at a temperature of about 1 K, and it does not 
seem to bear a direct relation to the critical temperature 
Tc 4 K) of the superconducting electrodes. We ob- 
serve that the crossover displayed in Fig. [1] corresponds 
to a temperature of about 1.5 K, in fair agreement with 
the position of the drop measured experimentally in the 
resistance. It is therefore quite likely that this feature 
may be an indirect signature of the temperature up to 
which the Cooper pairs are able to propagate without 
disruption, along the 2.5 /zm-long path of the sample re- 
ported in Ref. |l2. 



The fast decay of the supercurrent computed from ([7]) 
can be traced back to the vanishing density of states of 
graphene at the charge neutrality point. Actually, the 
form of the propagator in is a direct consequence 
of the conical dispersion around the Fermi points of 
graphene. One can therefore expect important changes 
in the Cooper-pair propagation upon doping the electron 
system. This can be investigated formally by introducing 
a finite chemical potential fj,, with the aim of shifting the 
Fermi level away from the charge neutrality point. At 



/i > 0, the Dirac propagator can be written as|2C 
G(°)(k,w) = {uj + VF<T-k)[ 



—llP- + v\\<P' — ie 



S(uj ~ vfIM) ^, I. iM 

+^^^^^^^^e{^i-VF\^)] (8) 

which is a convenient alternative form of expressing the 
electron propagator 



G(°)(k,c.) = 



(£(k)-^) , e{ii-em 



uj - e(k) 



oj - e(k) 



(9) 



for a dispersion with two branches £(k) — ±Wi?|k|. The 
Cooper-pair propagator can be computed now from the 
convolution of two Dirac propagators like ([8]). The result 

at a; = is 
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1 



_, , k arcsm — 

/i if |k| > 2fx 



(10) 



The slight change produced by the chemical potential 
in the infrared behavior of I?*^°^(k, 0) is enough to mod- 
ify the long-distance decay of the supercurrent. We have 
represented in Fig. [2] the result of evaluating 



from Eq. ([7|) with the Cooper-pair propagator in ([TT 
It can be seen that, at L ^ vp/lJ', the power-law decay 
of the critical current changes from to The 

supercurrents cannot be enhanced anyhow up to the mag- 
nitudes that they reach in a ID electron system, where 
the decay is given by a 1/L dependence in the noninter- 
acting theory. However, it may be worthwhile to explore 
experimentally the consequences of doping the graphene 
layer, up to levels where it can be affordable to measure 
supercurrents of the order of ~ 1 nA for suitably small 
values of L, according to the plot shown in Fig. [2] 

We also remark that doping the graphene layer may be 
the way to obtain experimental signatures of dynamical 
superconducting correlations. The pertinent approach to 
address this question is to deal again with the theory at 
nonvanishing chemical potential, starting from the Dirac 
propagator in Eq. ([8|). Computing now the Cooper- 
pair propagator at real frequency uj, the possible super- 
conducting instabilities have to reflect as singularities in 
D{0,uj), as it happens with the formation of bound elec- 
tron pairs in the BCS theory [ilj. 

The evaluation of the Cooper-pair propagator at chem- 
ical potential ^ > gives the result 
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FIG. 2: Plot of the zero-temperature critical current jP^^(O) 
(in units of lO'^evpkc ~ 1.2 ^A) as a function of the distance 
L (keeping W = lOV^c ~ 50 nm), for chemical potential 
11 ^ 1 meV (full line) and fi = 0.01 eV (dotted line). The 
dashed line corresponds to a decay and it is drawn here 
as a reference of that power-law behavior. 



We may sum up again multiple scattering processes of the 
Cooper pairs in the framework of the T-matrix approxi- 
mation. We observe then that D{0,uj) has a pole under 
conditions where the dominant interaction is attractive, 
such that V < 0. The denominator 1 -I- VD'-"'^ {0,u!) in 
([S]) vanishes at values of the frequency ui given by the 
equation 



|2/i — u!\ = u! exp ' 



ujV 



(12) 



There is always a solution of Eq. for uj below and 

close to 2/i when V < 0. Such a frequency marks the 
formation of bound electron pairs, that may take place 
for arbitrarily small values of /i, despite the fact that the 
density of states of graphene vanishes at cj = 0. The 
main difference with a conventional Fermi liquid is that 
the energy scale in front of the exponential in Eq. IT^ 
may be rather small, of the order of ^ 2fi. 

In conclusion, we have studied the proximity effect in 
graphene when the distance between the superconduct- 
ing contacts is much larger than their width. We have 
seen that, in the case of undoped graphene, the supercur- 
rents have a fast spatial decay, proportional to W^/L^. 
This strong dependence on L reminds of the behavior of 
the critical current in long diffusive junctions [22]. On the 
other hand, it is unlikely that the experimental measures 
may be strongly affected by disorder in a graphene layer, 
since backscattering arises only from small-range scat- 
terers with a size not larger than the lattice constant [isl]. 
We have actually seen that the results reported in Ref. 
12l are consistent with a crossover to the strong decay of 
the critical current at a temperature T* vp/kBL, and 



not at the much smaller scale given in a diffusive junction 
in terms of the diffusion constant P by T* ~ V/kBL^- 

We have seen that shifting the Fermi level away from 
the charge neutrality point changes the dependence of 
the critical current into a behavior, opening the 

possibility to observe a supercurrent over length scales 
above 1 //m at suitable doping levels. We have shown 
that this becomes feasible anyhow below the crossover 
temperature T* marking the onset of the power-law de- 
cay of the supercurrent, and that corresponds to the scale 
above which the Cooper pairs are increasingly disrupted 
by thermal effects during their propagation. 
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